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ABSTRACT 

We have developed the first self-consistent model for the observed large-amplitude oscillations along filament 
axes that explains the restoring force and damping mechanism. We have investigated the oscillations of multiple 
threads formed in long, dipped flux tubes through the thermal nonequilibrium process, and found that the oscillation 
properties predicted by our simulations agree with the observed behavior. We then constructed a model for the 
large-amplitude longitudinal oscillations that demonstrates that the restoring force is the projected gravity in the 
tube where the threads oscillate. Although the period is independent of the tube length and the constantly growing 
mass, the motions are strongly damped by the steady accretion of mass onto the threads by thermal nonequilibrium. 
The observations and our model suggest that a nearby impulsive event drives the existing prominence threads 
along their supporting tubes, away from the heating deposition site, without destroying them. The subsequent 
oscillations occur because the displaced threads reside in magnetic concavities with large radii of curvature. Our 
model yields a powerful seismological method for constraining the coronal magnetic field and radius of curvature 
of dips. Furthermore, these results indicate that the magnetic structure is most consistent with the sheared-arcade 
model for filament channels. 
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1. INTRODUCTION 

High-cadence Ha observations of large-amplitude standing 
oscillations along a filament axis were first reported by Jing 
et al. (2003), who measured bulk displacements of up to 
140 Mm, with a period of 80 minutes and a high velocity 
amplitude of 92kms _1 . Prominence oscillations are classified 
as “large amplitude” when velocity is of the order of 20 km s' 1 
or larger (Oliver & Ballester 2002). Similar large-amplitude 
longitudinal (LAL) oscillations were observed later by Jing 
et al. (2006) and Vrsnak et al. (2007). All the LAL oscillations 
reported so far show periods of 0.83-2.66 hr, damping times of 
2-10 hr, damping time per period of 2. 3-6. 2, and amplitudes of 
30-100 km s . These oscillations were apparently triggered by 
an energetic event: a sub-flare, a microflare, or a flare close to 
the filament. In addition, prominences undergo large-amplitude 
transverse oscillations in response to waves emanating from 
distant flares (e.g., Gilbert et al. 2008; Tripathi et al. 2009). The 
work presented here addresses only the longitudinal motions. 

Several models have been proposed to explain the LAL oscil- 
lations. Jing et al. (2006) interpreted their observed oscillations 
in terms of the Kleczek & Kuperus (1969) model. In this model, 
the prominence is considered a slab with a magnetic field along 
the filament, the restoring force is the magnetic tension, and the 
motion is perpendicular to the main axis of the filament. How- 
ever, this geometry is inconsistent with the observed direction of 
the LAL motion and the most likely magnetic structure: dipped 
or twisted field nearly parallel to the polarity inversion line (PIL; 
Mackay et al. 2010). The observations show that the displace- 
ments of the filamentary structures comprising the prominences 
are mostly aligned with the PIL and hence along the filament 
axes, without visible changes of the magnetic structure after 
the onset of the oscillation. In contrast, the magnetic tension 
drives only transverse motions, which would be perpendicular 
to the LAL oscillations. Vrsnak et al. (2007) assumed a flux- 
rope geometry and suggested that these oscillations are driven 
as a longitudinal-mode standing wave on a tightly coiled spring 


with fixed ends. This model also predicts motions perpendicu- 
lar to the local magnetic field, inconsistent with the observed 
LAL dynamics. Another possibility is that the restoring force is 
associated with gas pressure differences. However, the tempera- 
ture differences needed to accelerate the threads to the observed 
speeds are on the order of several million Kelvins, which are 
not observed (Vrsnak et al. 2007). The damping mechanism 
is also poorly understood. Several damping mechanisms have 
been suggested but not rigorously tested, e.g., energy leakage 
by the emission of sound waves (Kleczek & Kuperus 1969) or 
some form of dissipation (Tripathi et al. 2009; Oliver 2009). 

In this work, we propose explanations for the forcing and 
damping of the LAL oscillations based on our previous investi- 
gations of prominence formation and evolution (e.g., Antiochos 
et al. 2000). We found that all of the resulting condensations 
oscillate after their formation, and that these motions damp af- 
ter several oscillations. By studying in detail the motion and 
damping of many simulated threads in a sheared-arcade fila- 
ment channel (Luna et al. 2012), we could identify and analyze 
quantitatively the underlying physical processes. Our promi- 
nence model thus provides the first self-consistent explanation 
for the LAL oscillations. 


2. NUMERICAL MODEL AND RESULTS 

We focus our attention on the 47 individual oscillating threads 
modeled in our recent work (for details see Luna et al. 2012), 
which form within dipped flux tubes of a sheared-arcade fila- 
ment channel (DeVore et al. 2005); all flux tubes are assumed 
to remain rigid. The one-dimensional equations governing the 
plasma dynamics and energetics in these tubes are numerically 
solved with our ARGOS code (Antiochos et al. 1999). In re- 
sponse to steady, localized heating imposed in the tube foot- 
points, condensations of cool, dense plasma are produced in the 
corona by the well-investigated thermal nonequilibrium process 
(Antiochos & Klimchuk 1991; Antiochos et al. 1999, 2000; 
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Figure 1. Plot of the temporal evolution of the bulk velocity of the center of 
mass of the condensation example (solid line). The decaying exponential fit for 
the strong damping (dashed line) and weak damping (dot-dashed line) is also 
plotted (see the text). Additionally, the theoretical velocity fit obtained with our 
pendulum model described by Equation (5) is shown (dotted line). 



R(Mm) 

Figure 3. Scatter plot of the oscillation period of the threads as a function of 
the mean radius of curvature of the dipped part of the tubes (filled circles). The 
solid line shows the theoretical period from the pendulum model described by 
Equation (4). 



Figure 2. Scatter plot of the velocity amplitude as a function of the displacement 
amplitude for the set of 47 oscillating threads (filled circles), and observational 
data (diamonds) from Jing et al. (2003, 2006) and Vrsnak et al. (2007). 



Figure 4. Scatter plot of the strong damping time r s as a function of the weak 
damping time r w . 


Karpen et al. 2001, 2003, 2005, 2006; Karpen & Antiochos 
2008; Xia et al. 2011). 

We found that the combination of heating asymmetry and tube 
geometry determines the (usually off-center) position where the 
condensation forms, which is always closer to the footpoint with 
weaker heating, and the initial velocity of the thread. Threads 
rapidly settle down after oscillating around the bottom of the dip 
(Antiochos et al. 2000; Karpen et al. 2001, 2003), and accrete 
mass thereafter at a constant rate due to steady footpoint heating. 
As discussed in detail below, we find that the oscillatory features 
of the model prominence threads are strikingly similar to those 
of the LAL oscillations. The threads move mainly along the 
filament axis because the local magnetic field is highly sheared 
(i.e., nearly aligned with the PIL). 

A typical example of an oscillating thread is shown in 
Figure 1. Around t c = 2.7 hr, the thread forms with an initial 
velocity amplitude of 49 km s -1 but rapidly reaches a maximum 
speed of 53 km s _I . The maximum displacement in this case is 
40 Mm with respect to the bottom of the dip (the approximate 
equilibrium position). We note that the oscillation is very regular 
in time with a constant period of 1.0 hr. The oscillation is 
damped, but the velocity profile cannot be fitted with a sinusoidal 
function with a single decaying exponential envelope throughout 
the entire simulation. The initial oscillations damp rapidly at 
first followed by much weaker attenuation. We have fitted an 


exponential function for the first three periods of the oscillation 
envelope (dashed line), yielding strong damping with a decay 
time of r s = 2.3 hr. We have also fitted a decaying exponential 
function for the envelope of the curve after t = 10 hr (dot- 
dashed line), yielding a weak damping time of r w = 13.1 hr. 
The existence of these two damping times indicates the existence 
of two physically distinct phases in the temporal evolution of 
the thread. 

Similarly, we have computed the oscillation parameters 
for the 47 threads identified in our simulations. Figure 2 
indicates that a higher displacement amplitude is associated 
with a higher thread velocity. Most of the thread oscillations 
have velocity amplitudes larger than 20 km s -1 , and hence are 
classified as “large-amplitude” oscillations. We have not found 
any dependence of the periods on the length of the flux tubes 
or the mass of the condensations. However, we have found a 
dependence on the average radius of curvature of the dip, R, as 
shown in Figure 3. In Section 3 we explain the reason for such 
a dependence. There are few cases that do not fit very well the 
model curve. These tubes have a irregular shape and the radius of 
curvature is not uniform in the dip. The strong damping time r s 
as a function of the weak damping time r w is plotted in Figure 4. 
Clearly r s < r w even in the cases where the weak damping time 
is small. The decay time per period is r S /P = 1-6 for the strong 
and r w / P — 2-27 for the weak. 
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The clear agreement between the observed range of ampli- 
tudes, periods, and damping times (Section 1) and the results 
of our simulations (see Figure 2) led us to speculate that the re- 
ported LAL oscillations could be described with our theoretical 
model. We propose the following scenario to explain the ob- 
served periodic displacements, quasi-aligned with the filament 
axis, of parts of the filament. The LAL oscillations are triggered 
by an energetic perturbation at one end of the filament (e.g., Jing 
et al. 2006); the nature of the trigger is discussed in Section 5. 
This perturbation pushes the threads to oscillate at the same 
time, so the LAL oscillations represent the collective motion of 
bundles of threads that initially move in phase. In our model, the 
threads oscillate along the dipped part of magnetic flux tubes 
that are more or less aligned with the axis of the filament. Only 
a few LAL cycles were observed after the trigger, which we 
associate with the strong initial damping of the oscillations of 
our simulated threads. The weak damping has an small contri- 
bution to the attenuation in the first cycles. In addition, after 
several oscillations the excited threads lose their coherence due 
to slight differences in the oscillation periods, explaining why 
different parts of the filament oscillate increasingly out of phase 
as observed by Vrsnak et al. (2007). 

3. OSCIL L ATION PHYSICS: RESTORING 
FORCE AND DAMPING 


the motion is well described by the center-of-mass velocity. 
Although the threads periodically expand and contract slightly 
throughout the simulations, these perturbations do not affect 
the motion of the center of mass. We also assume that the 
amplitudes of the thread displacements are small with respect to 
the curvature radii of the dips, consistent with our simulations. 
The condensation exchanges mass and momentum with the 
surrounding corona, thus introducing the effects of drag into 
the force balance. Thus, we obtain the equation of motion as 


d 2 s 

dt 2 


+ 


1 ds 
t -t' dt 


+ or s = 0, 


( 2 ) 


and 


, m 0 tAo 

t — to = t 0 , 

a a) 


( 3 ) 


where s(t) is the position of the center of mass of the condensa- 
tion along the tube, a> is the angular frequency of the oscillation, 
and the phase i jf 0 = co m o /a . Equation (2) is similar to the equa- 
tion of a simple gravity-driven pendulum, with an extra drag 
term (second term). The coefficient of this drag term is the mass 
growth rate per unit of mass, ( dm/dt)/(m{t )). For the initial 
stages of the evolution ( t > t') the drag is intense, but for times 
much larger than t' this term is negligible. 

The angular frequency of the oscillation is 


We noted the existence of thread oscillations in previous nu- 
merical studies (e.g., Antiochos et al. 2000; Karpen et al. 2001, 
2003), but did not explore the origin of the restoring force and 
damping. Two possible origins for the restoring force have been 
proposed: pressure gradients and gravity. For several reasons, we 
find that pressure forces cannot account for the observed or sim- 
ulated behavior. Small-amplitude prominence oscillations were 
studied in two-dimensional magnetic configurations with slight 
flux-tube curvatures by Oliver et al. (1992, 1993), Joarder & 
Roberts (1993), and Oliver & Ballester (1995), who concluded 
that the oscillations were pressure driven. However, it is unclear 
whether these linear studies are applicable to large-amplitude, 
low-frequency oscillations. In the pressure-driven case, rapid 
changes in thread mass and length would cause frequency shifts 
(e.g., Morton & Erdelyi 2009; Erdelyi et al. 2011) that are not 
observed in either the solar or the simulated LAL oscillations. 
Finally, our simulations reveal that the pressure force on the 
condensation is negligible in comparison with the gravity force. 
Therefore, we conclude that pressure-driven oscillations are in- 
consistent with both the LAL observations and with the behavior 
of our simulated threads. 

We propose instead that the restoring force is mainly the solar 
gravity, and the damping mechanism is mainly the ongoing 
accretion of mass by the condensation. In the dipped part 
of a flux tube, the projected solar gravity is always directed 
toward the bottom of the dip. Therefore, the thread oscillation 
resembles the motion of pendulum with a length equivalent to 
the average radius of curvature of the dip, R. The temporal 
evolution of the cool mass (see Figure 3(a) of Luna et al. 2012) 
consists of an initial transient phase of strong growth starting 
at time t c and a linear regime after to- The transient phase is 
short (to — t c < P/2), so the mass of the condensation is well 
described as 

m(t) = m 0 + a(t - t 0 ), (1) 

where a is the growth rate of the condensation mass and mo is 
the mass at the beginning of the linear regime at to. In addition, 
all parts of the condensation move with similar velocities and 



where P is the period and go is the gravitational acceleration 
at the solar surface. Equation (2) is the zero-order Bessel 
equation with solution Jo[co(t — to) + j/rol- This oscillating 
function describes very well the first stages of the evolution: 
the amplitude decreases with time but not exponentially, which 
explains why we cannot fit a simple decaying exponential to the 
velocity profile shown in Section 2. After several periods the 
Bessel function is basically a sinusoid. However, the simulations 
also exhibit much slower damping during the later stages of the 
oscillation (see Figure 1), so it is necessary to introduce an 
additional damping term. The thread displacement takes the 
form 

s(0 = 50 + 71 Jo[mo(t - to) + fo)e- ( '- ,o),T - . (5) 

The velocity derived from this equation is overplotted in 
Figure 1 (dotted line). We obtain a very good match between 
the simulated oscillation and the oscillation model, which is 
also valid for the oscillations of the other 46 simulated threads. 
The additional weak damping might be associated with non- 
adiabatic effects such as conduction and radiation (Terradas 
et al. 2001; Soler et al. 2009). These losses are very effective 
in straight flux tubes where the curvature and gravity are not 
included, and the restoring force is the pressure gradient. In 
those systems, the non-adiabatic losses weaken the restoring 
force in each oscillation, damping the movement of the thread. 
However, in our system the main restoring force is the gravity 
projected along the tube, and the non-adiabatic losses are not so 
effective. 

Equation (5) indicates that both mass accretion and non- 
adiabatic losses contribute to the strong damping at the begin- 
ning of the oscillation. However, the former is more important 
because the latter term has a long damping time. The effective 
decay time associated with the mass loading is implicit in the 
Bessel function because its amplitude decreases rapidly in the 
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first few cycles. This damping is related to the phase \[r 0 of 
Equation (5). For i/c 0 = 0 the damping is maximum and the 
damping time takes its shortest value, but for larger values of 
\jfQ the damping is weak and the Bessel function is basically 
a sinusoid. Because this phase is inversely proportional to the 
mass accretion rate a (see Equation (3)), rapid mass loading im- 
plies strong damping of the oscillation and vice versa. However, 
relatively large values of the loading rate imply that i^o 0, 
indicating that the effective damping time is not very sensitive 
to a variations. The damping depends also on mo, which is not 
straightforward to characterize with a simple analytic approach. 
For these reasons, it is difficult to find a simple relation between 
the damping time and the parameters of the tube and the heating. 
This relation will be the subject of a future study. 

4. PROMINENCE SEISMOLOGY 

In our prominence model the magnetic structure exists in- 
dependent of the condensations. The shear and bending of the 
coronal magnetic field are determined by the shear at the pho- 
tosphere and the tension of the overlying field along the PIL 
(DeVore & Antiochos 2000). The magnetic-field structure is 
self supporting, and the interaction with the plasma content is 
small. However, the magnetic structure must support the denser 
threads. Thus, the magnetic tension in the dipped part of the 
tubes must be larger than the weight of the threads or equiva- 
lently 

B 2 

— - m n go > 0, (6) 

K 

where B is the magnetic field strength at the bottom of the dip, 
n is the particle number density in the thread, the mean particle 
mass m — 1 .27 m p (Aschwanden 2004), and m p is the proton 
mass. Combining Equations (4) and (6) we obtain 


combined with the large observed thread displacements, argue 
strongly for the applicability of the sheared-arcade magnetic 
structure. Assuming that the density of the observed prominence 
threads is in the range of the values found in our simulations, 
we estimate the minimum field strength in the prominence as 
approximately 34-68 G and 23-46 G, respectively, for the two 
cases. 

5. THE OSCILLATION TRIGGER 

In our simulations, steady footpoint heating produces con- 
densations that oscillate after formation. Each condensation has 
an initial velocity imparted by the thermal nonequilibrium pro- 
cess and determined by the heating asymmetry and flux tube 
geometry. However, the observed LAL oscillations clearly are 
related to an energetic event (subflare, microflare, or flare) near 
the filament ends. Although a detailed study of the excitation 
mechanism is beyond the scope of this Letter, we propose that 
the energetic event increases the heating impulsively at one end 
of the filament, increasing the evaporation of chromospheric 
plasma into the affected flux tubes. Our preliminary studies 
have revealed that condensations are remarkably robust; in- 
creased heating at one footpoint produces a flow of hot evap- 
orated plasma that pushes a thread to a new location, rather 
than destroying it in situ. For example, in the Jing et al. (2003) 
event the subflare pushed an entire section of the filament 70 Mm 
away from the initial position. We speculate that extra heat is de- 
posited only on those footpoints located closest to the energetic 
event, so only parts of the filament are excited simultaneously 
and oscillate collectively in phase. In addition, the increased 
evaporation produces more mass accretion and increases the 
oscillation damping. 

6. CONCLUSIONS 


B ^ 


go m » „ 

4t r 2 ’ 


(7) 


which constrains the minimum field strength as a function of 
the thread oscillation period. Equation (7) can be written as 

B[GI;S2<i (wT^) ,,!,>|hrl - <8) 

For the thread electron number densities of our simulations 
in = 1.6-6. 2 x 10 11 cm -3 ) and the oscillation periods ( P = 
0.75-1.75 hr). Equation (8) yields a minimum field strength 
of B Js 31— 75 G for the 47 threads with a mean value of 
5 1 G, in agreement with the sparse observations of intermediate- 
type prominence magnetic fields (Mackay et al. 2010). In our 
prominence model, the deep-dip flux tubes have an average area- 
expansion factor of three with respect to their footpoints. Thus, 
we estimate the mean field strength at the flux-tube footpoints 
as 153 G, also consistent with observations (e.g., Aschwanden 
2004). 

In addition to providing valuable information on the difficult- 
to-measure coronal magnetic field in filament channels, our 
model can also determine the radius of curvature of the promi- 
nence flux-tube dips from the observed oscillation properties. 
We have used the Jing et al. (2003) and Vrsnak et al. (2007) ob- 
servations to solve Equation (4) for R. The observed periods are 
80 minutes and 50 minutes, respectively, yielding correspond- 
ing curvature radii of 152 Mm and 62 Mm. This large-scale 
curvature of the field deduced from the observed oscillations, 


We conclude that the observed LAL motions represent a 
collective oscillation of a large number of cool, dense threads 
moving along the magnetic field, triggered by a nearby energetic 
event (subflare, microflare, or flare). The main restoring force is 
the projected gravity in the flux tube dips where the threads 
oscillate, even in threads with large radii of curvature. The 
period is independent of the tube length and of the constantly 
growing mass. The oscillations are strongly damped by the mass 
accretion of the threads. After several oscillations the excited 
threads lose their coherence due to slight differences in the 
oscillation periods, causing different parts of the filament to 
oscillate increasingly out of phase as observed. The mechanism 
that initiates the collective oscillation is not determined in the 
present investigation. We speculate that the energetic event 
impulsively increases the evaporation at the filament flux- 
tube footpoints closest to the energetic event, increasing the 
evaporation rate only on one side of the threads. The hot 
evaporated plasma pushes the existing threads, initiating the 
oscillations. A detailed study of the initiation mechanism will 
be the subject of future work. More observations of oscillating 
prominences are needed to further refine and test our model. 
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